Analysis of a SU(4) generalization of Halperin's wave functions as an approach 
towards a SU(4) fractional quantum Hall effect in graphene sheets 
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Inspired by the four-fold spin- valley symmetry of relativistic electrons in graphene, we investigate 
a possible SU(4) fractional quantum Hall effect, which may also arise in bilayer semiconductor 
quantum Hall systems with small Zeeman gap. SU(4) generalizations of Halperin's wave functions 
[Helv. Phys. Acta 56, 75 (1983)], which may break differently the original SU(4) symmetry, are 
studied analytically and compared, at v = 2/3, to exact-diagonalization studies. 
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The fractional quantum Hall effect (FQHE) in conven- 
tional semiconductor (typically GaAs) heterostructures 
may be understood to great extent with the help of one- 
component U(l) trial wave functions, such as Laughlin's 
[l| or composite- fermion generalizations of it c 2:]. The 
effect occurs due to the repulsive electronic interactions 
when the filling factor, the ratio v — n e i/riB between the 
electronic density n e i and the flux density ub = B/(h/e), 
belongs to the "magic" series v = p/(2sp± 1), with inte- 
gral s and p, or particle-hole symmetric fillings. 

Soon after Laughlin's original proposal, Halperin 
pointed out that the one-component picture for the 
FQHE may have some shortcomings due to a relatively 
small Zeeman effect, Az ~ 0.3_B[T]K, which docs not 
outcast in GaAs the characteristic interaction energy 
e 2 /cIb ~ 50y^3[T]K until extremely high magnetic fields 
Instead, he proposed a two-component generaliza- 
tion of Laughlin's wave function in order to account for 
the spin SU(2) degree of freedom Indeed, polar- 

ization measurements indicate that several states, such 
as e.g. v = 2/3 and 2/5, are not spin-polarized and 
may sometimes undergo transitions between competing 
FQHE states with different spin polarization when the 
total magnetic field is varied at constant filling factor 0] . 

The recent observation of a relativistic integral quan- 
tum Hall (QH) effect in graphene @, Q has led to several 
theoretical investigations concerning an eventual FQHE 

in this novel carbon compound m B EE EEI Q. 

In contrast to the abovementioned GaAs heterostruc- 
tures, graphene has an underlying SU(4) symmetry due 
to the two-fold valley degeneracy in addition to the 
physical spin degree of freedom. The Coulomb interac- 
tion respects this symmetry to lowest order in a/Is ~ 
10 _2 ...10 _1 , where a = 0.14nm is the distance between 
nearest- neighbor carbon atoms, and Is = y/h/eB = 
25/y A B[Tjnm is the magnetic length [81. fill]. In previous 
exact-diagonalization studies on potential candidates of a 
FQHE in graphene UGH, the SU(4) symmetry has been 
omitted and a complete spin polarization presupposed. 
In this case, the two-fold valley degeneracy may be mim- 



icked by a SU(2) isospin, and graphene may be treated 
as a non-relativistic SU(2) QH system if one replaces 
the non-relativistic by the relativistic effective potentials 
[H, [H[ ■ More recently, Toke and Jain have proposed a 
composite-fermion construction at v = p/(2sp + 1) with 
an internal SU(4) symmetry 

Here, we investigate a graphene SU(4) FQHE, where 
we omit the assumption of complete spin-polarization 
because the ratio A z /(e 2 /el B ) ~ 0.003. ..0.014-/B[T]K 
remains small (the precise value depends on the effec- 
tive dielectric constant which may vary from e ~ 1 for 
a free graphene layer to e ~ 5 on a SiO substrate 
We propose a SU(4) generalization of Halperin's wave 
function, which may lead to a plethora of new FQHE 
states not captured in previous studies. We mainly dis- 
cuss the filling factor v = 2/3, where our SU(4) exact- 
diagonalization results indicate rich physics beyond the 
SU(2) case 0, Notice that the graphene filling fac- 
tor vq is defined with respect to the Dirac points, and 
the central Landau level (LL) is therefore half-filled at 
vq = 0, whereas in semiconductor QH systems v is de- 
fined with respect to the bottom of the lowest LL. In 
order to make the connection between the two of them, 
taking into account the four-fold degeneracy, one has to 
choose vq = —2 + v (or vq = 2 — v, due to electron-hole 
symmetry) . 

In the spirit of Ref . Q , we generalize Laughlin's wave 
function [lj to SU(-ftT), with K components, 
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TABLE I: Examples of SU(4) wave functions. 



takes into account correlations between the different com- 
ponents. Here, denotes the complex position of the 
fcj-th particle of component i, the exponents rrij and 
indicate the strength of the intra- and intercomponent 
(i =/= j) correlations, respectively. Similar generalizations 
have been suggested by Qiu et al. for multilayer QH sys- 
tems [l5[ and by Morf in the framework of the FQHE 
hierarchy scheme We concentrate on K = 4, with 

1 = (T, +), 2 = (t, -), 3 - (I, +), and 4 = (|, -), where 
T, l denote the ^-component of the physical spin and ± 
the valley in graphcnc or the layer index in GaAs bi- 
layer systems. We (artificially) break the spin-valley [2(| 
symmetry and distinguish between intracomponent pseu- 
dopotentials 17 1 V^, V^~ s for different spin in the same 
valley, and V£ for electrons in different valleys. Notice 
that the wave function |T]) is an exact ground state for 
a model interaction > for m < min(mi, ...,rri4,), 
V r f t > for m < min(rii2, ni4, n 2 3, n 34 ), V^~ s > for 



m < min(ni3, 77,24), and V r 



A/E/E-s 



otherwise. In 



GaAs bilaycr QH systems, the difference between inter- 
and intralayer interactions fixes V£ = Vj^~ s . 

For general K, the component filling factors Vj = 
Pj/n>B of j-particles may be obtained from the usual zero- 
counting argument according to which the flux density 
n b equals the sum of the j-component densities pj times 
the strength of their zero, fi^ = Pj m j + Pi n iji f° r 

all j. With the help of the symmetric K x K exponent 
matrix Mk = { n ij), where riji = riij and rijj = rrij, one 
may thus calculate the v.j in a concise manner, 



if the matrix Mx is invertiblc 



Mk 1 (1, ■ 



1 
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Tab. I shows examples of SU(4) A^-particle 
wave functions ([1]), labeled by the set of exponents 
[mim2?7i3m4, n e n+n_], where we have restricted the 10 
exponents to only 7, n + = 7713, n_ = 7724, and n e = 
7712 = flu = ^23 = ^34- This means that we treat all 




FIG. 1: The three polarizations, S Z ,I Z , and P z , are restricted to 
the interior of a tretrahedron, completely explored in the case of 
a fully SU(4)-symmetric state (r = 1). For states with a matrix 
rank r = 2, the polarizations are restricted to quadratic planes with 
fixed S Z ,I Z or P z , e.g. [3333,033] state with I z = (plane in light 
gray). The [3333, 111] state, with r = 4, has S Z ,I Z ,P Z = (black 
dot). 



intervalley correlations on the same footing, with an ex- 
ponent n e , and the intravalley correlations for different 
spin are described by the exponent n+ and ?i_ in the val- 
leys + and — , respectively. Apart from the total filling 
factor vt = v\ + V2 + ^3 + v&, the [7771777277737774, 77 e Tj+77_] 
wave functions have the spin S z j (N/2) = [v\ +i>2 — (vz + 
Vi)]/vT-, the valley I z /(N/2) = [vi +1/3 - (v 2 + v A )]/v T , 
and a third polarization P z /(N/2) = \v\ + v± — (i>2 + 
v, i)\l v T- These polarizations arc good quantum numbers 
associated with the mutually commuting z-components 
of the SU(4) spin, t z ®1, 1®t z , and t z ®t Zi respectively, 
in terms of the diagonal Pauli matrix t z . Although they 
vary in the interval [— N/2, N/2], they are not indepen- 
dent and indeed restricted to the interior of a tetrahedron 
depicted in Fig. [T] 

Some physical properties of the SU(4) states may be 
characterized by the rank r of the matrix M4. In the 
case of r = 4, as e.g. for the [3333, 111] wave function at 
v = 2/3 or [5555,222] at v — 4/11, all Vj are uniquely 
determined by Eq. ([2]), and the S Z ,I Z , and P z are thus 
fixed. This is represented by the black dot in Fig. Q] 
for the [3333,111] wave function. For r = 3, three rows 
of M4 are linearly independent. Consider the first row 
[(Ti+) component] to be a multiple of the third [(!,+) 
component]. In this example, the wave function, such as 
[3737,233] at v = 4/11, represents a state with a SU(2) 
spin ferromagnet in the + valley Defining a combined 
filling factor for the + component, v + = + , one 
may describe this state alternatively by a SU(3) wave 
function {1}, with an invertable M3 matrix. For r = 2, 
one may e.g. realize a SU(3) ferromagnet of three com- 
ponents, whereas the combined three-component filling 
factor and that of the fourth are fixed. Another possi- 
bility is that the SU(2) ferromagnet discussed for r = 3 
is accompanied by another SU(2) ferromagnet in the — 
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valley with no coherence between the two of them. Exam- 
ples are the [3333, 033] and the [3535, 235] wave functions 
at v = 2/3 and 4/11, respectively. The polarizations are 
now restricted to planes with fixed I z , as depicted in Fig. 
[TJ By a simple exchange of the components, harmless 
in the case of an underlying SU(4) symmetry, S z or P z 
may play the role of I z so that the planes with fixed 
I Z ,S Z , and P z are equivalent. Both ferromagnetic states 
for r = 2 may be described by a SU(2) Halperin wave 
function with invertable matrix M2. The case r = 1 rep- 
resents a SU(4) ferromagnet the polarizations of which 
explore the full tetrahedron depicted in Fig. [T] The only 
constraint is fixed by ut, and one finds a SU(4) Laughlin 
wave function, where all possible states may be obtained 
by any SU(4) rotation of a state with V\ = l/(2s + 1) 
and V2 = V3 = ^4 = 0. 

As becomes apparent from Tab. I, the same filling 
factor may be realized by different wave functions with 
different matrix rank r. This has to be contrasted to 
the SU(2) case where one needs to invoke the composite- 
fermion theory Q to obtain competing states with dif- 
ferent spin polarization at the same filling factor, such 
as e.g. at v = 2/5. However, filling factors which do 
not arise in U(l) or SU(2) wave functions, such as e.g. 
v = 8/19 (or 10/23 and 26/47, not shown in Tab. I), 
may be described by SU(4) wave functions with r = 4. 
One may therefore principally expect a closer vicinity of 
FQHE states in SU(4) than in SU(A" < 4). 

In our numerical studies, we use an exact diagonal- 
ization on the sphere with a fully implemented SU(4) 
invariance. All calculations are performed in the lowest 
LL. The various SU(4) trial wave functions are obtained 
by tuning the pseudopotentials in the abovementioned 
manner. They appear as zero-energy ground states with 
the lowest number of flux quanta for a given interaction 
(similarly we obtain the Laughlin state or the Halperin 
states |18j|). The number of flux quanta threading the 
sphere is 25. Due to the large Hilbert spaces we need 
to consider, our calculations are restricted to the system 
sizes of N — 8 fermions. 

In the spherical geometry, the two states [3333,111] 
and [3333,033] occur at 25 = (3N/2) - 3. Fig. E^a) 
shows the energy spectrum of the FQHE at v = 2/3, 
with 25 = 9 and N = 8, for a choice of pseudopoten- 
tials V* = V® = V£- s = 0, except for Vf = V E = 
Vq~ s = 1, for which [3333,111] is expected to be the 
exact ground state. The ground state is found at a to- 
tal angular momentum L = 0, which is, together with 
the finite gap to all excited states, a condition for an in- 
compressible FQHE state. Moreover, the ground state is 
non-degenerate and has S z = I z = P z = 0, in agreement 
with the [3333, 111] state (Fig. Q). In Fig. HJb), the en- 
ergy spectrum (25 = 9, N = 8) is shown for a different 
pseudopotential choice with V£ = = Vj%~ 3 = 0, ex- 
cept for Vq~ s = Vf~ s = = 1 - a model for which 
the [3333, 033] state is expected to be the exact ground 
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FIG. 2: Low-energy spectra, for N = 8 and 25 = 9, with contact 
potential in the LLL, built in such a way that the trial states (a) 
[3333, 111] and (b) [3333, 033] are the only ground states with zero 
energy. We only show the lowest-energy state for each value of the 
total angular momentum L. 
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FIG. 3: Low-energy spectra with Coulomb interaction in the lowest 
LL (a) for 25 = 9 and N = 8; (b) for 2S = 11 and N = 8. We 
show only the lowest energy for each value of L in (a), but also 
the spectrum for 5 Z = 0, I z = 0, P z = 4 (crosses), for comparison 
with the SU(2) case [sj. The continuous line shows a 1? fit through 
the lowest-energy states. The only fit parameter is known with a 
precision of 1% (if removing the L = 1 point) in both cases. 



state. One finds a gapped L — ground state with I z = 
and a degeneracy due to a free choice of S z and P z within 
the light gray plane in Fig. [1] as for a [3333,033] state. 

The Coulomb interaction in the lowest (central) LL 
turns out to be more involved. We have first performed 
Monte-Carlo calculations for up to N = 70 particles, 
which indicate that in this case the [3333,111] state is 
much lower in energy than the [3333, 033] state when 
compared to the statistical error for both states. Fig. 
[3£a) shows the energy spectrum for 25 = 9 and N = 8, 
obtained from exact diagonalization of the lowest-LL 
Coulomb interaction. The spectrum is reminiscent of 
that for 25 = 11 and N = 8 [Fig. Mb)], but it is clearly 
different from those of the [3333,111] and [3333,033] 
states. Both spectra show a L — ground state that 
is connected to a collective mode with a striking L 2 be- 
havior (black lines). However the ground-state degener- 
acy is different, as shown in Tab. II. For 25 = 11, we 
obtain a [N/2, N/2] = [4,4] multiplet [H], in agreement 
with results by Toke and Jain [12j. The non-degenerate 
S z = I z = 0, P z — 4 state, is, up to an interchange of S z 
and P z , the state described by Apalkov and Chakraborty 
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TABLE II: Degeneracy of the ground state for the Coulomb inter- 
action in the different < S z < I z < P z sectors for N = 8, 25 = 11 
(second row) and N = 8,25 = 9 (third row). (004) and (013) are 
not ground states for 25 = 9. 



who have considered an internal SU(2) valley symme- 
try and a complete spin polarization. Indeed, the spectra 
coincide, but Fig. [3jb) shows clearly that this is only the 
top of the iceberg and that states with different SU(4) po- 
larization also belong to the ground-state manifold. More 
importantly, the low-energy properties are dominated by 
a L 2 mode (Fig. [3]) we conjecture to be associated with 
collective spinlikc waves the symmetry of which is deter- 
mined by a subgroup of SU(4). Notice that such mode 
would not affect the incompressibility of the ground state. 



In order to make a connection between the [3333, 111] 
ground state and that obtained for Coulomb interaction 
in the lowest LL for 2S = 9 and N = 8, we evalu- 
ate the ground states when varying the pseudopotcn- 
tials V E and from (exact model for the [3333, 111] 
state) to 1, keeping Vi — V E = 1. Here, we set 
V E = V^~ s for all m. Fig. [5] shows the resulting 
phase diagram, in terms of the ground state degenera- 
cies. The [3333,111] state is stable over a large region 
in the phase diagram, but one also obtains doubly de- 
generate (at moderate V® ~ 0.4... 0.5) and compressible 
ground states with L ^ 0. Most saliently, a 105 x degen- 
erate ground state is obtained at V® = 1, for V% > 0.25, 
such as that corresponding to the Coulomb interaction 
(dashed line in Fig. [4]). This degeneracy is precisely that 
of a [4,4] SU(4) multiplet. The phase is critical in the 
sense that it is destroyed as soon as we slightly deviate 
from V E = Vf = V* = 1. 

In conclusion, we have analyzed a SU(4) generalization 
of Halperin's wave function [3j, which may be a promiss- 
ing approach towards a SU(4) FQHE eventually occuring 
in high- mobility graphene sheets. These trial wave func- 
tions yield incompressible states at filling factors, which 
may not, in some instances, be described by Laughlin 
[H, SU(2) Halperin [jj, or SU(4) composite-fermion wave 
functions [13], such as at v = 8/19. Moreover, the rank 
of the exponent matrix may be related to ferromagnetic 
properties of the trial states within the SU(4) symmetry. 
Whether these FQHE states are realized, depends on the 
precise form of the interaction potential. At v = 2/3, 
two candidates, [3333, 111] and [3333, 033], with different 
residual symmetry, yield incompressible FQHE states for 
appropriately chosen pseudopotentials. The Coulomb in- 
teraction, however, chooses a state with a pronounced 
L 2 collective (Goldstone) mode, both for 2S — 9 and 11. 
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FIG. 4: Left panel : schematic view of the ground state degeneracy 
for pseudopotentials VP and VP varying from to 1, evaluated on 
a regular grid (11 X 11 values), for 25 = 9,7V = 8. The white 
region corresponds to a non-degenerate ground state, light grey 
regions are doubly degenerate ones with S z = I z = P z = and 
dark grey regions are 105 X degenerate ground states. Hatched 
regions are compressible ground states with L ^ 0. The black dot 
corresponds to the ratios VP = Vf /V^ and Vp = V-f /V^ , where 
C 

V m ' are Coulomb pseudopotentials in the lowest LL. Right panel : 
gap as a function of VP for fixed VP = 0.4. The peak corresponds 
to the 105 X degenerate ground state. For Vp > 1, there is a 6 
times degenerate ground state with one of the quantum number 
maximaly polarized and the others unpolarized. 



This mode is thus likely to survive in the thermodynamic 
limit. 
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